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Inelastic neutron scattering from superconducting (SC) Ba(Feo.926Coo.o74)2As2 reveals anisotropic 
and quasi- two- dimensional (2D) magnetic excitations close to Qafm = (i|) - ^he 2D antiferromag- 
netic (AFM) wave- vector of the parent BaFe2As2 compound. The correlation length anisotropy of 
these low energy fluctuations is consistent with spin nematic correlations in the J1-J2 model with 
J1/J2 ~ 1- The spin resonance at ^8.3 meV in the SC state displays the same anisotropy The 
anisotropic fluctuations experimentally evolve into two distinct maxima only along the direction 
transverse to Qafm above ^80 meV indicating unusual quasi-propagating excitations. 

PACS numbers: 74.20.Mn, 74.25.Ha, 74.70.-b, 78.70.Nx 



The unconventional superconductivity in the newly 
discovered i?FeAs(Oi_xFx) {R = rare earth) super- 
conductors [1] with critical temperatures up to 55 K 
[21 emerges upon suppression of an antiferromagnetic 
(AFM) phase and is accompanied by a spin resonance in 
the magnetic excitation spectrum [4-9 , indicating a close 
connection between magnetism and superconductivity. 

In the tetragonal structure of parent ferropnictides, 
degenerate magnetic states and frustration arise from 
competition between nearest-neighbor (NN) and next- 
NN (NNN) AFM exchange couplings. Magnetic frustra- 
tion is thought to produce an emergent nematic degree of 
freedom [10 , which couples to orthorhombic distortions, 
inducing a structural transition pTHT3] . Nematic fluctu- 
ations may also be present in the superconducting (SC) 
phase [13 , potentially leading to local spin and electronic 
anisotropics within the Fe layers. The relationship be- 
tween nematicity and unconventional superconductivity 
in iron arsenides [I4 is a subject under intense debate. 

Here we examine the in-plane wave- vector dependence 
of magnetic fluctuations in SC Ba(Feo.926Coo.o74)2As2 
[15] over a wide range of energies using inelastic neutron 
scattering (INS). Our results show that the anisotropic 
spin fluctuations within the Fe layer consistent with the 
square symmetry of the Fe layer evolve into unusual 
modes propagating only along the direction of ferromag- 
netic (FM) near neighbor spin correlations. We also find 
an analogous Q-space anisotropy of the spin resonance 
in the SC state. 

High-quality Ba(Feo.926Coo.o74)2As2 single crystals 
display a sharp SC transition at Tc ^ 22.2 K [15 . 
INS measurements were performed on the wide Angular- 
Range Chopper Spectrometer (ARCS) at the Spallation 
Neutron Source with incident energies up to 250 meV 



and on the HB-3 triple-axis spectrometer at the High 
Flux Isotope Reactor with fixed final energy at 14.7 meV. 
The mosaic of the coaligned samples (~5.40 g, 26 crys- 
tals) for the ARCS (HHL) measurements is ~3.30° and 
'^3.50° full width at half maximum for rotations about 
the (HHO) and (OOL) directions, respectively. The HB-3 
measurements in the {HKO) scattering plane were car- 
ried out on a realigned subset of six of the ARCS crys- 
tals with a total mass of ~2.58 g and a mosaic spread of 
-0.53° for both (i^OO) and (Oi^O) directions. 

An overview of the Q dependence of the magnetic scat- 
tering at 4 K measured by the ARCS is shown in Figs. 
[TJb) andjljc). Due to the 2D nature of the spin fluctua- 
tions [4, 5 , 16 , we present L-integrated data. Maxima are 
observed for (i^, K) = (m+^,n+^), m and n being in- 
tegers, as for the AFM stripe ordering in BaFe2As2. The 
scattering peaks are quite broad and feature a distinct 
anisotropy, indicating short-ranged spin correlations with 
an anisotropic correlation area in the paramagnetic state. 
Similar observations in Ba(Feo.935Coo.o65)2As2 [16^ and 
paramagnetic CaFe2As2 [17 suggest that this anisotropy 
is a universal property of spin fluctuations in the iron 
arsenides. The spin-space anisotropy of magnetic fluctu- 
ations in BaFe1.9Nio.1As2 [IE] bears not directly related 
to the reciprocal space anisotropy observed here. 

In Figs. [2|c)-[2|f), all scattering patterns display a 
twofold symmetry with respect to Qafm- particular the 
Q-width is considerably larger along the direction trans- 
verse (TR) to Qafm than along the longitudinal optical 
(LO) direction. Two distinct maxima split off in the TR 
direction at —100 meV as shown in Fig.[2|f). Figure [TJa) 
shows that the TR direction at (^^) corresponds to FM 
spin correlations in the stripe AFM structure. At first 
sight, it might seem that this twofold symmetry breaks 
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FIG. 1: (color online) (a) A single Fe layer showing a local spin 
configuration corresponding to Q^fm ~ (||) (NNN) 
exchange interactions Ji (J2)- The tetragonal (down) and Fe 
lattice unit cells are shown as dashed squares. Overview of 
the Qab = Q — (Q -0)0 dependence at 4 K measured on ARCS 
with (b) 50 meV and /kj = 20 ± 4 meV and (c) 

250 meV and /kj = 75 di 15 meV. An empty sample- holder 
background was subtracted. The dashed lines in [(b) and (c)] 
show the Brillouin-zone boundary of the Fe square lattice. 
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FIG. 2: (color online) Energy and wave- vector dependences 
around Qa6 = (||) [(a) and (c)-(f)] in the SC state and (b) 
in the normal state, (a) and (b) show the Q = (HKO) depen- 
dence at the resonance energy ^8 meV with the same inten- 
sity scale from HB-3. [(c)- (f)] Integrated data over (c) 17-24 
meV, (d) 40-50 meV, (e) 70-80 meV, and (f) 100-136 meV 
with (c) Ei = 50 meV and [(d)- (f)] 250 meV averaged over 
two quadrants from ARCS. An empty sample- holder back- 
ground was subtracted. 



the fourfold symmetry of the Fe square lattice. However, 
Fig. 1 shows that the (|^) point already has twofold 
symmetry. While the twofold symmetry of the observed 
scattering pattern thus does not break the symmetry of 
the Fe sublattice, it indicates a spin-correlation area that 
is anisotropic with respect to the direction of correlated 
NN spins. 

The quality of the data can be ascertained in constant 
energy cuts passing through (^^) along the LO and TR 
directions (circles in Fig. [3|. We analyze our data by 
starting with the assumption that, in the \ow-q^ low-o; 
limit, the spin fluctuations near an AFM critical point 
can be described by a diffusive model [19l [20^ . As in Ref. 
[17], this model can be extended to include anisotropic 
spin correlation lengths within the Fe layer, 

x''(Qafm + q.^) 



^2 + r2[i + + ^^e?R]^ ' 

Here xo represents the strength of the AFM response 
function, F is a damping constant and ^lo and ^tr are 
spin correlation lengths along LO and TR directions rel- 
ative to Qafm- I^i^s to the neutron intensities, S{Q^uj) = 
/'(IQI)x"(QAFM+q,w)(l-e-'^'^/''^)-^ where /2(|Q|) 
is the magnetic form factor of the Fe^+ ion, were per- 
formed by a procedure where ^ was varied while keeping 
F = 9.5 ± 1.0 meV as obtained from fitting the 30 K data 
in Fig. [4]^a). The resulting averaged anisotropy of the 
LO and TR spin-correlation lengths is given by the ratio 
^Lo(10.4 dz 0.6 A)/^tr(5.9 ± 0.4 A) = 1.8 ± 0.2 at 4 K. 
Despite the expectation that the diffusive model is valid 



only at low energies, we find that the parameters above 
describe constant energy cuts over the entire energy range 
in the LO direction, as shown in Figs. 3(a)-3(e). Along 
the TR direction, the diffusive model works well up to 
-80 meV [Figs. 3{d)-3{e')]. Above -80 meV, the TR 
spectrum splits into two peaks as shown in Figs. 3 (a') 
and 3(b') that cannot be described by the diffusive model 
(to be discussed). 

The spectral features within the diffusive model are 
described by a single peak centered at Q^pj^ whose q- 
space half-width at half maximum (HWHM) depends on 
energy according to <7hwhm — ^~^(V2 + CJ^F"^ ~ !)• To 
confirm the validity of the diffusive model over a wide 
energy range, we show that the Lorentzian HWHM of the 
constant energy cuts vs energy agrees with the ^hwrm as 
shown in Fig. 3(f), which also highlights the anisotropy 
between LO and TR widths. 

The in-plane anisotropy is consistent with an itinerant 
description of the dynamical magnetic susceptibility due 
to the ellipticity of the electron pockets at Q^pM [21] • Al- 
ternatively, a local-moment picture yields the same phe- 
nomenology and allows the anisotropy to be associated 
with competing NN (Ji) and NNN (J2) exchange interac- 
tions. We thus note that our analysis cannot address the 
question of the appropriateness of an itinerant or local- 
moment description for the iron arsenides. For a region of 
stripe AFM correlations corresponding to Qafm = (^|)' 
the LO (TR) direction is along AFM (FM) correlated 
spins [Fig. 1(a)]. When evaluated in the J1-J2 model, 
the dynamical magnetic susceptibility near 

^^AFM takes an 

identical form to model Eq. (1), allowing the anisotropy 
observed in the correlation lengths to be associated with 
the exchange ratio. Specifically it can be shown that 
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FIG. 3: (color online) [(a)-(e)] LO and [{a)-{e)] TR scans 
(circles) around Qab = (||) with Ei = 250 meV and energy 
transfers as indicated at 4 K from ARCS. The approximate 
L (r.l.u.) ranges are [(a) and (a )] [6.4, 7.5], [(b) and (b')] 
[4.6, 5.8], [(c) and (c^] [3.2, 3.8], [(d) and (d^] [2.7, 3.2], and 
[(e) and (e')] [2.2, 2.7]. The solid lines are fits of [(a)-(e) and 
(c^)-(e^)] the diffusive model and [(a^) and (b^)] the ballistic 
model. The dashed lines are calculated with [(a^) and (b^)] the 
diffusive model and [(c^)-(e^)] the ballistic model with fitted 
parameters held fixed. The horizontal bars in (a), (a^), (e), 
and (e^) represent the expected q resolutions, (f) compares 
the width of single Lorentzian fits along LO (circles) and TR 
(triangles) directions to the HWHM of the diffusive model 
(lines), (f) compares the position of TR Lorentzian peaks 
(squares) to the dispersion (slanted dashed line) and peak 
splitting of the ballistic model (solid line). The horizontal 
dashed line corresponds to the energy c/^e as described in 
the text. 



= 2 i°j| [TT. ^TR < ^Lo indicates that Ji is AFM and 

"'2 SlQ+StR 

acts to destabilize (stabilize) FM (AFM) near-neighbor 
correlations. Using the experimental correlation lengths 
we obtain ^ = 1.0±0.2at4K, which is similar to the ra- 
tio found in Ba(Feo.935Coo. 065)2 As2 [TF, and in the AFM 
ordered CaFe2As2 [22l[23], and indicates NNN interac- 
tions are important in the iron arsenides. We emphasize 
that the observed anisotropy does not break the four- 
fold symmetry of the Fe square lattice but it does imply 
anisotropic correlations and hence interactions between 
NN spins. The inferred exchange ratio is indeed within 
the regime of frustrated magnetism ( < 2) where spin 
nematic properties have been predicted. 

The absence of any splitting in the LO direction is 
unusual because one expects the diffusive excitations 
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FIG. 4: (color online) (a) Energy dependence of INS at Q = 
(||0) below (5 K) and above (30 K) Tc (solid symbols) after 
background subtraction. The solid lines are fits as described 
in the text, (b) Raw scattering data as well as their respective 
backgrounds taken at Q = (0.683 0.183 0) from HB-3. 



to evolve into damped spin-wave modes at short wave- 
lengths (high energies). Assuming a typical conical spin- 
wave dispersion, constant energy cuts should display 
ringlike features at high energy, instead of split maxima. 
Rather, the TR splitting bears some similarity to mag- 
netic excitations in incommensurate spin-density-wave 
systems such as Cr [24], and in the FeTei-a^Sea^ pnictide 
superconductors [2F. The phenomenological Sato-Maki 
function [26j has been used to describe the spin fluctua- 
tions in these incommensurate systems. However, in our 
case, there is no indication of incommensurability at low 
energies and we thus treat the TR excitations as damped 
propagating modes by including a ballistic term in the 
diffusive model. 



;2+r2(l+g2^2_^)2' 



(2) 



where c corresponds to the velocity of the propagating 
mode and is an effective TR correlation length. The 
TR splitting above '^80 meV can be fit to Eq. (2), shown 
as solid lines in Figs. 3(a') and 3(b'). This produces 
= 7.4 ± 0.8 A and c = 245 zb 10 meVA. The obtained Ie 
is somewhat larger than ^tr = 5.9 ± 0.4 A from Eq. (1). 
Despite the unusual splitting, we find that c is compara- 
ble to the spin- wave velocity in the TR direction in AFM 
ordered CaFesAss (300-350 meVA) [22l [23l. Lester et 
a/. [16 report a similar TR velocity (230 meVA) in para- 
magnetic Ba(Feo.935Coo.065)2As2. 

To better determine the evolution of the TR splitting, 
we also fit the experimental TR spectra with a symmetric 
pair of Lorentzians, whose peak splitting is plotted in Fig. 
3(f). According to the ballistic model, a diffusive (single- 
peaked) response is obtained below cjd = c^^^ ^ 33 meV 
in the TR direction. Above this energy, the ballistic 
model predicts split peaks at g = ±c~^^s/uj'^ — c^^~^, 
which only approaches a damped simple harmonic os- 
cillator (DSHO) response at high energies (with splitting 
q = :^uj/c)^ as shown in Fig. 3(f). The actual obser- 
vation of a splitting depends on the statistical quality of 
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the data which, in our case, ahows a clear observation 
of the TR sphtting only above ^80 meV. However, the 
relatively low value of uj^ explains the agreement between 
the TR velocities obtained from our ballistic model and 
DSHO model of Lester et al. [16] . In the LO direction, 
the magnetic response appears diffusive at all energies 
(subject to the statistical quality of the data), appar- 
ently due to a much higher LO velocity. The report of 
a LO dispersion with Clq = 580 meVA based on DHSO 
model analysis is probably a lower bound. The finite cor- 
relation length of the system allows one to estimate that 
= Clo^lo^ ^ meV suggesting that most if not all of 
the LO data analyzed in Ref. [16 is within the diffusive 
limit. 

Despite the agreement of the diffusive/ballistic mod- 
els along the LO and TR directions, the nature of the 
excitations in the full [H^ K] plane is still unclear. The 
unidirectional nature of the split modes and the strong 
anisotropy cannot be easily modeled by Eq. (2), which 
highlights the anomalous nature of the high-energy mag- 
netic spectrum. 

We now discuss the effects of superconductivity on the 
anisotropic spin excitations, which were examined using 
the triple- axis spectrometer HB-3. Figure [4] shows the in- 
elastic intensity at constant Q = (^^0) as a function of 
energy transfer at 5 and 30 K. Spectral weight is pushed 
to higher energies as superconductivity develops. To ob- 
tain magnetic scattering only, we measured identical en- 
ergy scans at a point in reciprocal space Q = (0.683 0.183 
0) [Fig. |4jb)], where no appreciable temperature depen- 
dence is observed between 5 and 30 K. We subtracted 
their average from the data at Q = (^^0) to obtain 
the background-subtracted magnetic scattering spectra 
shown in Fig. [4]^ a). The spectra can be fit, respectively, to 
a single imaginary-pole-response function in the normal 
state (overdamped response) and a damped harmonic os- 
cillator in the SC state. The normal state relaxation rate 
at the critical wave- vector is found to be Fn = 9.5 ± 1.0 
meV. In the SC state the damping rate has decreased to 
Fs = 6.6 ±0.4 meV and the resonance energy is hVt ~ 8.3 
meV. 

To examine the spatial correlations associated with the 
spin resonance, we measured the momentum dependence 
at 8 meV and 5 K [Fig. [2ja)]. A quasielliptical feature 
around Q = (^^0) was observed similar to the data be- 
low ^80 meV. The anisotropic scattering extends beyond 
the instrumental resolution ellipsoid shown in Fig. |2|a). 
At 30 K [Fig. [2|b)], the normal-state spin fiuctuations 
display a similar momentum space anisotropy. The LO 
and TR cuts through the resonance can also be well fit 
to the diffusive model yielding peak positions (TR) and 
peak widths that are consistent with those data as shown 
in Figs. |3|f) andjsjf). We conclude that within errors 
there are no changes in the spatial correlations above and 
below Tq. It is predominantly the spectrum of magnetic 
excitations that is modified by superconductivity. 



To summarize, we have observed collective magnetic 
excitations in SC Ba(Feo.926Coo.o74)2As2 close to the 
Qafm wave- vector of the parent BaFe2As2 compound. At 
low energies, the excitations have a pronounced in-plane 
anisotropy that can be associated with frustrated versus 
satisfied NN interactions in the parameter regime of dy- 
namic nematic correlations. The spin resonance in the SC 
state is found to have the same anisotropy. At energies 
above ~80 meV, quasipropagating modes are observed 
experimentally along the TR direction, while modes in 
the LO direction appear to have a very large energy scale. 
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